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We discuss averaged turbulence modeling of multi-scales of length for an incompressible 
Newtonian fluid, with the help of the maximum information principle. We suppose that 
there exists a function basis to decompose the turbulent fluctuations in a flow of our 
concern into the components associated with various spatial scales and that there is a 
^ , probability density function / of these fluctuation components. The unbiased form for / is 

determined and the turbulence model is closed, with the multi-scale correlations up to the 
c/3 ■ fourth order, through maximizing the information under the constraints of equality and 

inequality for that flow. Due to the computational difficulty to maximize the information, 
a closely related but simple alternative objective is sought, like the determinant or the 
trace of the second order correlations of the turbulent flow. Some preliminary results and 
Oh implications from the application to homogeneous turbulence are presented. Some issues 

yet to be resolved are indicated. 

5^ 1 Introduction 

There are two widely used methodologies in Reynolds averaged turbulence modeling: One 
^ ! adopts multi-point multi-scale correlations of fluctuations, while limited to homogeneous tur- 
O ! bulence (cf. [I] and [ID]); The other, without the restriction of homogeneous turbulence, adopts 
one-point correlations and builds widely used engineering models (cf. [8] and [9]). There are 
also developments of multi-point or mixed closure models for inhomogeneous turbulence (cf. [2] 
^ . and pj]). Exploring a different strategy, Tao, et al. (cf. [12]) have attempted to unify the 
methodologies by both employing multi-scale correlations to inhomogeneous turbulence and 
resolving the issue of closure through multi-objective optimization under various constraints 
of equality and inequality. These constraints are either general and intrinsic to turbulence or 
special to a particular flow concerned. 

Specifically, Tao, et al. (cf. [12]) have presented a rationale for statistically averaged multi- 
scale turbulence modeling for an incompressible Newtonian fluid on the basis of optimizing 
certain objective functions motivated by homogeneous turbulence. From the viewpoint of 
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information availability, they have justified why their formulation is restricted to the resolution 
of the average flow fields up to the second order correlations and to the modeling of the third 
and fourth order correlations in algebraic forms. Also, they have explained why they choose an 
optimization approach to help make the average model determinate, which can be understood 
from the need to satisfy the constraints as to be listed in Section [2j 

Several objective functions are tentatively proposed to be minimized in [T2]. One difficulty 
encountered in the proposal is the development of a computational scheme to carry out the 
multi-objective optimization itself and the selection of an optimal solution among the possibly 
many; the alternative is to seek a single objective function to be optimized as discussed in 
this work. Another difficulty faced is that there appear no clear-cut rules to guide in the 
determination of the turbulence model construction. We attempt to resolve these issues in this 
work, with the help of the maximum information principle (cf. [6] and [7]). We should also 
point out that some objectives proposed in [12] are incorrect, such as the minimization of the 
fluctuation kinetic energy. 

One major issue in turbulence modeling is the lack of data on higher order correlations, e.g., 
about the boundary conditions for the one-point correlations higher than the Reynolds stress 
in an inhomogeneous turbulent flow. This is the issue of information unavailability. One of 
its implications is that we encounter difficulties to test, calibrate or apply a turbulence model 
containing higher order correlations as part of its primary fields. Also, a model containing 
higher order correlations faces the formidable challenge in the aspect of mathematical analysis 
and numerical computation due to the great number of variables, constraints and evolution 
equations involved. Therefore, we intend to construct a model including and resolving only 
the lower order correlations of the velocity and pressure fluctuations, and thus, the model 
cannot be expected to provide accurate data for the correlations beyond. If we construct 
such a model, it is informationally sound to have a corresponding probability density function 
which may produce adequate results for the included correlations with the maximum 'spread- 
out' allowed, i.e., the probability density function is unbiased. The maximum information 
principle offers such a tool, which is pursued in this work. Moreover, the principle apparently 
leads to the maximum entropy production in the aspect of the micro-states defined by the 
velocity fluctuations of various spatial scales. The principle has single objective function, the 
information, to be maximized, a preferable alternative to multi-objective optimizations. 

We are aware of the controversy and the problem surrounding the maximum information 
principle, especially in the case of continuous probability distributions. We should also mention 
the question of the applicability of the principle to turbulence model construction, considering 
the non-equilibrium nature of turbulent motion. Here, we intend to explore the possibility and 
consequence of applying the principle to turbulence modeling and to offer an interesting test of 
the principle itself. We will see that the adoption of the maximum information principle leads 
to a closed multi-scale turbulence model up to the fourth order correlations. There still remains 
a tremendous task to make the model computationally feasible to turbulence simulation, as to 
be delineated later. 

For an objective function involving the probability density function explicitly, such as the 
information in the maximum information principle, there appears an intrinsic difficulty to 
make the objective optimization computationally feasible which will become clear in Section [3j 
This motivates us to seek some closely related but computationally less demanding alternative 
objective functions to be maximized, like the determinant or the trace of the second order 
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correlations and so on. 

To check whether the present model may produce meaningful results, we apply it to the 
special case of homogeneous turbulence; some preliminary results are given so as to indicate its 
potential. 

The present work intends to explore a non-traditional methodology to resolve the issue of 
closure and construction of multi-scale turbulence models. It indicates, in details, the challenge 
to construct an averaged turbulence model resolving the multi-scale correlations of velocity 
fluctuations and pressure fluctuations, in comparison with engineering turbulence modeling. 

This paper is organized as follows. In Section [2J we present the basic ideas about the 
turbulence modeling of multi-scales and the constraints of equality and inequality for the fluc- 
tuation correlations up to the fourth order. We explore in Section [3] the maximum information 
principle and its consequence in turbulence modeling, with the information as the single ob- 
jective function to be maximized. In Section HJ as a simplification of and an approximation to 
the information, we discuss the invariants of the second order correlation as single objectives. 
In Section [51 we apply the model to homogeneous turbulence and discuss some preliminary 
results and implications. Some concluding remarks are made in Section [6] about the present 
formulation and some of the issues to be investigated. 

2 Basic Formulation 

Consider the isothermal motion of an incompressible Newtonian fluid in a finite domain 
T> C M 3 at time t > 0. The fields for the velocity and pressure are governed by the Navier- 
Stokes equations, i.e., 

vj,j=0, (2.1) 

dvi , . . . 

— + [ViVj)^ = -p H + . (2.2) 

Here v is the kinematic viscosity of the fluid; the mass density p is incorporated into p. Adopting 
the standard practice in Reynolds averaging, we take the partition 

Vi = Vi + Wi, p = P + q, w~i = Q, q = 0, (2.3) 

where Vi and pP are the ensemble averaged velocity and pressure which behave regularly; Wi 
and pq are the parts associated with the random fluctuations of the velocity and pressure fields, 
if the flow is turbulent. Then, (12. 1 p through (12.31) result in 

Vj, 3 = 0, (2.4) 

dVi , , 

-r£ + (ViV 3 + nW])v = -P,i + v V hjj , (2.5) 

w jv = 0, (2.6) 
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dt 
and 



dw ■ 

1 + (wiWj + WiVj + ViWj - Ww]),j = -q,i + v , (2.7) 



{wiWj + 2wiVj - WiWj),ij . {21 



Next, as done in [12], we assume that there is an orthonormal basis for the functions defined 
in£0 

B := |fr(x; m) : x G V, m = (mi, m 2 , m 3 ), m^ G {0} UN, J 6(x; m) 6(x; n) c/x = 5 m n |, 

(2.9) 

so that the fluctuations Wi and q can be adequately represented as 

Wi = Wi(t; m) 6(x; m) , g = q(t; m) 6(x; m). (2.10) 

m m 

Each m may be interpreted as a wave number characterizing certain spatial scales, along with 
its associated 6(x; m) which is oscillatory in space similar to that of Fourier representation; 
Wi(t; m) and pq(t; m) are the fluctuation components of velocity and pressure associated with 
the scale of m. From the consideration of the deterministic continuum treatment and the 
viscous dissipation effect, there is expectedly an upper bound on m, and thus, the number of 
the terms in the above representations is finite. It is easy to extend the function basis to the 
case in which a weight is required to define orthogonality. It follows from (12.61) through f l 2 . 1 f) 
that the evolutions of the second and third order correlations of 



Aj(iJ) :=*i(*;i)^'(*;j)> 7ijfc(i,j,k) := Wi(t; i) Wj(t; j) w k (t; k) (2.11) 
as well as the correlations of pressure fluctuations are governed by 

^Ai(U) - ^ £(m; i) /%(m, j) + ^ (w a (m, i) + W(ia; i) <*«)/?y(m, j) 

m m 

D(m; j) /^(m, i) + ^ (Wji(m, j) + W(m; j) fy) /3 K (m, i) 

m m 

A,(m; n, i) 7^(n, m, j) - £;(m; i) g(m)^-(j) 

m,n m 

^A,(m;n,j) 7^(11, m,i) - 5j(m; j) g(m)u>;(i), (2.12) 



In the case of T> being a cuboid, we may use the Cartesian coordinate system and the Legendre polynomials 
or the Chebyshev polynomials, etc. If the turbulence is homogeneous, one may resort conventionally to Fourier 
series. There are cases in which a curvilinear coordinate system and other function basis may be preferred. For 
a rather generally shaped 2?, we might still employ the Legendre polynomials or the Chebyshev polynomials 
as the function basis (in the Cartesian coordinate system with a cuboid containing £>), however, the functions 
defined on V need to be extended to the cuboid properly. 
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^2 D ( m 'i P) q(m)Wi(i) = 2 ^ Wfc( m ; P) /?fci( m > i) + C '^( m ' n ' P) 7fc*i(m, n, i) , (2.13) 



d_ 
dt 



7^(i,j,k) -u^D (in; i) 7 ijfc (m, j, k) + ^ (^( m > + W(m; i) <5 a W fc (m, j, k) 

m m 

i/^D(m;j)7jjfei(m,k,i) + ^ (w}i(m,j) + W(m;j) ^^(m, k, i) 

m m 

^£>(m; k) 7 fcij (m,i,j) + ^ fw w (m,k) + W(m; k^W^m, i,j) 



v 

m 



- - ^ Ai(m; n, i) ^^ fc (n, m, j, k) - /3 K (n, m) /3 jfc (j, k) j - ^ ^(m; i) g(m)tfy(j)iuj(.(k) 

m,n m 

- y^A;(m; n,j) f<$y W (n, m, k, i) - /%(n,m) /3 fei (k, i) j - Bj(m; j) g(m)«5 fc (k)«Ji(i) 

m,n m 

- y^A(m; n,k) (<J IWi (n,m, i, j) - Afc(n,m) /%(i,j) J - ^ -Bfc(m; k) g(m)wi(i) Wj(j) , 

m,n m 

(2.14) 

- D & i) g(l)^(j)w fc (k) = 2 ^ W,(m; i) 7 y*(m, j, k) 

1 m 

+ 5^C/i(m,n;i) foy fc (m, n, j, k) - ^«(m, n) /3 jfe (j, k)Y (2.15) 

m.n 

and 

- ^ £>(1; i) ?(l)g(p) = 2 £ W,(m; i) g(p)^(m) + ^ C, 7 (m, n; i) g(p)^(m)^(n) . (2.16) 

1 m m,n 

Here 



<Jy«(i, j, k, 1) := Wi(i) Wj(j)w k (k)wi(l) (2.17) 

is the fourth order correlation, and 

Wi(m) := Wi(t; m) , g(m) := q(t; m) , 6(m) := &(x; m) ; 

A, (m;n, k):= / [6(m)],j 6(n)6(k)dx, £>j(m;n): = / [6(m)],j 6(n)dx, 
Jv Jv 

C^(m,n;k): = / [6(m)]„- [6(n)],< 6(k)dx, D(m;k):= / [&(m)],ji 6(k)dx, 
H/(m;k):= / [b(m)],j Vj 6(k)dx, Wi(m;k):= / [6(m)]„- 6(k)dx, 
Wij-(m,k):= / 6(m)6(k)V r i , i dx. (2.18) 



The Reynolds stress tensor in (j2.5p can be represented as 



um- = J2&Ahi)Ki) &Ci)- (2.19) 

Equations (12.41) . (12. 5ft . ( 12.12ft through ( 12. 19ft form a closed model of multi-scale turbulence, 
provided that the fourth order correlations are approximated, like in terms of the lower order 
quantities. From the solutions of this model, we may construct the multi-point correlations 
such as, V x, y, z, a G V, 

viifatywjfy^) = ^/%(i,j)6(x; i)6(y;j), 

ij 

Wi(x, t) Wj(y, t) w fc (z, t) = ^2 7y fe (i, j, k) 6(x; i) 6(y; j) 6(z; k), 

ij,k 

u?i(x, t) iwj(y, t) io fe (z, t) «>,(a, t) = ^ S ijk i(i, j, k, 1) 6(x; i) 6(y; j) 6(z; k) 6(a; 1), 

ij,k,l 

?(x,t)g(y,t) = E^HO) 6(x; k) 6(y; 1). (2.20) 

k,I 

The inclusion of (I2.14p for the evolution of 7r,-fc(i, j, k) is based on the consideration that 
7jjjk(i, j, k) and 5ij k i(i, j, k, 1) are needed to obtain the correlation of pressure fluctuations as 
seen above. If we model j, k, 1), it is preferable to have an equation of evolution to 

describe 7ijjt(i, j, k); the treatment preserves the intrinsic relationship between 7^(1, j,k) and 
5ijki(i, j, k, 1) as displayed in ( 12. 14ft and (12.151) and the constraints to be discussed below, and 
also, it makes the fourth order correlation the sole quantities to be modeled. 

This discussion justifies the adoption of a multi-scale turbulence model up to the fourth order 
correlation. Models involving higher order correlations are not pursued due to the information 
unavailability and the difficulty of mathematical analysis and computation. These problems 
are encountered even in the present model, especially regarding the inclusion of (I2.14p : The 
first is how to furnish adequate initial conditions for 7^(1., j, k) in the case that the average flow 
fields are time-dependent - its resolution may be case dependent. The second is how to furnish 
boundary conditions for WjWjWk, in the case that we do not have such data - we leave their 
determination as part of the solution from the model through the optimization to be introduced 
below. 

The question now is how to determine <%/w(i, j, k, 1) such that the model composed of (12.41) . 
(12. 5p . ( I2.12p through (1 2 . 1 9 [) is closed. In the following, we explore specifically the scheme of the 
maximum information principle so as to construct a determinate model through maximizing 
one objective function called the information. To this end, we discuss first the constraints 
of equality and inequality on /3^-(i,j), 7yfe(i, j, k) and j, k, 1) as an integral part of the 

formulation of closure and maximization. 

i. The definitions of the correlations imply that 

ftj(iJ) = &i(i> i) > 7ijfc(U, k ) = 7 jifc (j, i, k) = j kji (k, j, i) = 

5 ijk i(i, j, k, 1) = 5 jik i(j, i, k, 1) = 8 kj u(\s, j, i, 1) = 5i jki {l, j, k, i) = ... (2.21) 
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That is, each correlation remains invariant under the interchange of any two pairs of its 
indexes, say {i, i} and as indicated above. 

ii. Applying the Cauchy-Schwarz inequality to 7yfc(i,j,k) and ^^(i, j, k, 1) gives, 
Wi, j, k, I, Vi, j, k, 1, (the summation rule suspended), 

i) > 0, i) j) - (/%(i, j)) 2 > 0, 

i, j, j) > 0, i) *jifcfc(j, j, k, k) - (y ijk (i, j, k)) 2 > 0, 

^(i,i,j,j)4^(k,k,l,l) - fe,(i,j,k,l)) 2 > 0. (2.22) 

iii. The Reynolds stress tensor should be positive semi-definite, i.e., 

WiWi > 0, WiWi w 2 w 2 — {wiw-z) 2 > 0, det [wjifj] > 0, (2.23) 

in order to guarantee the non-negativity of both the fluctuation energy and the average 
bound for the velocity fluctuations. 

The quantity of Wi,k Wj,k should be positive semi-definite, 

wi, k wi,k > 0, w 1 , k wi,k w 2 ,iw 2 ,i - (w 1 , k w 2 ,k) 2 > 0, det[wi,kWj,k\ > 0, (2.24) 

due to its significance associated with the viscous dissipation in turbulence. In symmetry, 
one may also enforce the positive semi-definiteness of Wk,i Wk,j similarly. 

Considering the importance of vorticity in turbulence (cf. [21), one may introduce Uij : = 
(wi,j —Wj,i )/2 and formulate the constraints of 

(^F>0, (^P>0, JuuF>0, (^(^-(cJT^) 2 >0, 

{lj 23 ) 2 (o; 3 i) 2 - (u 23 u 31 ) 2 > 0, {u 31 ) 2 (u 12 ) 2 - (uj 31 u 12 ) 2 > 0, (2.25) 

from the definition and the Cauchy-Schwarz inequality. 
It is also natural to require that 

g2>0, q 2 wf - (qmtf > 0, q 2 {wiWj) 2 - {quWjf > 0- (2.26) 

iv. Equations (12. 6p and (12. 10 j) result in, Wi, j, k, Vi, j, k, n, 

^ Am(i, m) B m (m; n) = 0, ^ Jijm{hi, m) B m (m; n) = 0, 

m m 

Sijkmih j, k, m) B m (m; n) = 0. (2.27) 

m 

v. We apply the Cauchy-Schwarz inequality to the multi-point correlations of fluctuations 
in the physical space so as to obtain the following constraints, i.e., Wi, j, k, I, m, n, p, 
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Vx, y, z, a G V, (the summation rule suspended), 



(wi(x)) 2 (wj(y)f - (wi(-x.)wj(yj) > 0, (c^(x)) 2 (u H (y)) 2 - (u)ij{x) w«(y)) > 0, 
(u7 i;m (x)) 2 (wj, n (y)) 2 - (wi, m (-x)wj, n (y)\ > 0, 

2 2" 2 /" \ ^ 

(lUj(y)tUfc(z)) > 0, (Wi(x)) (u7j(y) w; fc (z)) J - ^ i (x)w i (y)w; A .(z)J > 0, 

(io»(x) Wj(y)) 2 (iUfc(z) w;( a )) 2 - (™;(x) Wj(y) Wfc(z) iui(a)) > °> ( 2 - 28 ) 

and so on, including more spatial derivatives of first order or even higher orders. One 
may also obtain the constraints involving the pressure fluctuations. 

vi. There are constraints from the boundary conditions for the Reynolds stress tensor on dV, 
whose forms and number are flow-dependent; there is also the possibility of zero WjWjuJk 
and WiWjWkWi on the solid wall boundary of a flow if the wall is smooth and experiences 
negligible random vibrations. 

If one has additional information about a specific flow, one may formulate it as additional 
constraints holding in that flow. 

Several comments are in order regarding the above constraints, (i) The constraints formu- 
lated in the physical space reflect partially the link between the equations of evolution ( 12. 12ft 
through (12.161) and the counterpart equations of evolution for it?j(x, t) Wj(y, t) and so on in the 
physical space, (ii) The constraints of inequality from the Cauchy-Schwarz inequality are not 
that tight since the inequality holds for any two integrable functions. As a part of the model, the 
collective effect of the very many such constraints may be sufficiently severe to the correlations, 
(iii) The constraints of ( I2.23P through (I2.26P and ( I2.28P hold in T> and the constraints from the 
boundary conditions hold on &D, they need to be treated appropriately in order to derive the 
corresponding constraints for /^-(ij), 7yfc(i, j,k) and <%jy(i, j, k, 1). One possible scheme is to 
integrate the inequalities (holding in T>) on the domains of T>, DxD, etc. with weight functions 
like (6(x;m)) 2 , (6(x;m)) 2 (&(y;n)) 2 , etc. (iv) The number of the constraints is great, and 
thus, the approximate treatment of representative wave numbers should be adopted to make 
the model computationally feasible, as to be discussed in Section [3J Also, the approximate 
treatments of locality and local isotropy to be discussed will convert part of equations (I2.12p 
through (12.161) into constraints of equality, (v) There is also the issue about Kolmogorov's 
results as well as their modifications (cf. [1] and [5]). We need to examine whether these results 
are the natural outcome of the model with the constraints listed above or the results themselves 
should be explicitly formulated and imposed as additional information and constraints in the 
model. 

3 Maximum Information Principle 

To apply the maximum information principle (MIP), we introduce a vector w of components 
{wi(t; k)} which are the turbulent velocity fluctuation components of all the spatial scales, at 
instant of t, involved as in (I2.10p . We also introduce a probability density function / of the 
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velocity fluctuation components, /(w). The introduction of this probability density function 
is justifiable in that we have effectively used it in setting 

/ /(w)dw = l, / i2J i (i)/(w)dw = tw i (i) = 0, / Wj(j) f(w)dw = j) , 

JH»N J^SN J R 3N 

/ wfi) Wjd) w k (k) f(w)dw = 7y fc (i, j, k) , 

JR3N 

/ Wi(i) Wj(j) u)fc(k) u)j(l) /(w)dw = <J ijW (i,j,k,l) , (3.1) 

where A" is the number of the wave numbers {m} involved and 3A" the dimension of w. 

Additional remarks are worth making about the form of or the variables involved in / (w) . 
Clearly, this is a partial treatment, in contrast to a formulation of the probability density 
functional for Vi and p (cf. [1]) and from it the derivations of (12. 4p . (12. 5p . (I2.12p . (I2.13p . etc. Here, 
we do not pursue such a mathematically comprehensive and rigorous treatment. We suppose 
that the turbulent velocity fluctuations can be defined as in (12. 3p from the notion of ensemble 
average, these fluctuations can be further decomposed into the components associated with 
various scales of length according to (12.101) . and a probabilistic description of these fluctuation 
components can be carried out. The pressure fluctuations q(m) can be represented in terms of 
tfr as seen from ([21]), fl2~T0l . (T2~T3|) . (12~T5D and f l2~T6|) . The effect o f the average velocity and 
pressure on / will be accounted for through certain parameters or coefficients contained in /, 
as to become clear below. 

We may formulate the principle in the form of 



/ f($)\nf($)d* (3.2) 



where I is the information to be maximized (cf. [B] and [7]). Since w may be viewed as 
representing micro-states of turbulent fluctuations, the maximization of / is like the maximum 
entropy production regarding the micro-states, if all the constraints on the flow are known and 
imposed. 

Now, we maximize (13. 2p under (13. ip . with the help of the method of Lagrange multipliers, 
to yield 



/(w) = exp 



A - ^2 M*) ^i(i) _ M'(i> j) ^i(i) ~~ M'fe(i' k ) w<(i) ^/(j)™fc( k ) 

(3.3) 



- ^2 ^ijki(i,J,Kl) Wj(j)wk(k)wi(l) 

ij,k,l 

Here, A, Aj(i)'s, A^- (i, j)'s, A^^i, j, k)'s and \ijki(i, j, k, l)'s are the Lagrange multipliers enforc- 
ing, respectively, the five sets of constraints in (13. ip . 

All the fourth order correlations are yet to be fixed and expected to be determined with 
the help of the maximum information principle through maximizing /. Substituting (13. 3p into 
(13. 2p . we have 

I = A + ^2 M(i> j) fti(iJ) + M fe (i> J' k ) 7ijfc(iJ> k ) + A ^(i> J, k > 4/fc/(i, j, k , !)■ 

ij ij,k ij,k,l 

(3.4) 
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The multipliers of A's and the correlations up to the fourth order are linked through (13. ip . 
Theoretically, we may search and determine the values of Xijkiih j, k, l)'s by maximizing /. 
Specifically, we find these values, for a specific flow of our interest, by maximizing I under 
(12. 4p . (12.51) . (I2.12p through (I2.19p and (13.11) . the aforementioned constraints of equality and 
inequality, the initial conditions for the second and third order correlations, and the initial and 
boundary conditions for the average velocity as well as the boundary condition for the average 
pressure. That is, we may solve the turbulent flow of interest. 

We notice that the number of the correlations, the multipliers, the equations and the con- 
straints is of astronomical order of magnitude, and in practice, we have to reduce it significantly 
with the help of certain approximations physically plausible. Also, the maximum information 
principle is conceptually based on inference in the face of information incompleteness, and more 
data from the physical consideration will help to narrow the spread-out of / and to improve the 
turbulence modeling. There are several conventional approximations we may adopt to these 
ends: 

(a) Instead of employing a complete function basis, we will thin the function basis through 
selecting representative wave numbers and the associated basis functions sparsely and 
appropriately, and thus, we can reduce the number of the correlations as well as the mul- 
tipliers to be solved; for instance, we may take only the wave numbers m = (m 1 ,m 2 , ^3) 
whose components m 1; m 2 and m 3 correspond to nearly the same physical length scales. 
This treatment had better be implemented in the representations of equation (I2.10p so as 
to derive the equations of evolution for /^-(ij) and 7yfc(i,j,k) and the equations for the 
correlations of pressure fluctuations conveniently and consistently. 

(b) We resort to the standard notion of locality or cascade such that all the correlations 
between the wave numbers not neighboring to each other are neglected. This point needs 
to be clarified further by quantifying the meaning of neighborhoods. 

(c) To adopt the conventional treatment of small-scale, locally isotropic turbulence, we define 
the notion of high wave numbers and treat the velocity fluctuation components associated 
with these high wave numbers as isotropic. For example, we set /3 ij7 - ( i , j ) = j) 5^, if 
the associated wave numbers are high. 

The approximations above offer some possibilities to reduce the complexity of the closed 
model, but they do not provide great help to the computation of (13.11) which is formidable due 
to the high dimensions involved. It is desirable to have a simplified alternative version of (13. 2p 
which does not involve / explicitly, as to be discussed next. 

4 Approximate Objective Function 

It is known that the information I for a Gaussian is proportional to the determinant of the 
covariance matrix j)] (cf. [5]); Also, from probability theory, the covariance matrix is a 

major indicator regarding the extent of the spread-out of /, 




(4.1) 
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The mathematical properties of [/^(i, j) ] are characterized by its non-negative eigenvalues 
and equivalently by its invariants; Two invariants of interests are the trace ^ k (3 kk (\i, k) and the 
determinant det j) ] . In the case that the probability density distribution of a turbulent 
flow is slightly non-Gaussian, the above-mentioned result of [3] implies that 

7 D :=det[/%(i,j)] (4.2) 

is an adequate objective function, alternative to the information /. The difficulty of such 
a choice is the relevant computational complication. From the perspective of computational 
convenience, the use of 

/T:=$^Mk,k) (4.3) 

k 

as an objective is appealing. Another appeal of this choice is its physical significance, consid- 
ering that the total fluctuation kinetic energy possessed by a turbulent flow in T> is related to 
It through 

If 1 . 1 

- / «»rfx=-V/3 tt (k,k) = -/ T (4.4) 
Z J ^ i i 

(in the case that the function basis is orthonormal) . The third appeal is the relatively simple 
structure of the trace in the space of {5ij k i(i, j, k, 1)} defined by 

5 ijkl (i,l k, 1) := 6 ijk i(i, j, k, 1) - /%(i, j) /3 k i(k, 1) 

which is motivated by the structures of (I2.14p and f)2.15p . 

From the perspective that turbulence arises from instabilities due to the change of some flow 
parameters in flows which are originally laminar, such as the increase of Reynolds number, a flow 
tends to have non-trivial turbulent fluctuation energy under proper flow conditions (cf. [13]). 
This does not necessarily imply that It of (14. 3[) is the maximum for the flow to occur under 
the constraints and conditions given; this point may be clear if one takes, say Id of f ]4.2[) as 
the objective function to be maximized. Further analysis and specific simulations have to be 
done in order to decide which of the invariants is suitable as the alternative objective function 
approximating the information I. 



5 Homogeneous Turbulence 

To test preliminarily whether the turbulence model constructed above has the potential to 
produce meaningful results, we apply it to the special case of homogeneous turbulence here. 
Let us consider a homogeneous turbulent motion in which Vi = Vy Xj in the cube of T> = [L/2, 
—L/2] 3 with Vij being constant and V kk = 0. In addition to the constraints of equality and 
inequality listed in Section [2J there are constraints of equality such as 

(wi(x) Wj(x)), m = 0, (wi, k (x) Wj,l(x)\ m = 0, (wi(x)wj(x) io fe (x)J, m = 0, 

(wi(x.) Wj(x) w k (x) wi(-x)\ m = 0, (5.1) 
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that are special to homogeneous turbulence resulting from its supposed spatial homogeneity. 
We may integrate these constraints of equality in the domain of T> with weight functions like 
6(x; m) in order to derive the corresponding constraints of equality for j), 7yfc(i, j, k) and 
Sijki(i, j, k, 1). As a consequence of (I5.ip . the WjWj related terms disappear from (12.51) . (12. 7p 
and (12. 8p . and the product /3y(i, j) /3jw(k, 1) related terms disappear from (I2.14p and (I2.15P ; the 
average flow field of V{ and P is not affected by the correlations. To discuss the consequence of 
these disappeared terms clearly, we gather all the relevant evolution equations here: 

^Aj(U) -v^D{m-\) /%(m, j) + (w«(m, i) + W(m; i) <J«) /3y (m, j) 

m m 

■ ]T £>(m; j) /^(m, i) + £ fe(m, j) + W(m; j) fy) /3„(m, i) 



m 



m,n m 

- J^i4 i (m;ii,j)7y i (n,m,i) - ^^(m; j) q(m)wi(i), (5.2) 

m,n m 

- ^ D ( m 'i P) ?( m )wi(i) = 2 Wfe( m ; P) /9w(m, i) + C< ^( m ' n 5 P) 7fcZi( m > n, i) , (5.3) 

m m m.n 

^7y fc (i,j,k) - i/^D(m; i) 7y fe (m, j, k) + ^ (Wfl(m, i) + W(m; i) 7^(111, j, k) 

m m 

i/^D(m;j)7jjfei(m,k, i) + ^ (w#( m J) + ^( m ; j) *ji)7iw(m, k, i) 

m m 

uy^D(m; k)7 Wj -(m,i,j) + ^ (w k i(m,k) + W(m; k) t^W^m, i, j) 

m m 

= ~ y^A t (m; n, i) <% fc (n, m,j,k) - ^^(m; i) g(m)to i (j)t(} fe (k) 

m.n m 

- A/(m; n, j) 6i jki {n, m, k, i) - ^ £j(m; j) g(m)u)j(.(k)u)i(i) 

m.n m 

- J^Aj(m;n,k) 5/ fc y(n, m, i, j) - ^ -B fc (m; k) g(m)tBj(i) Wj(j) , (5.4) 

m.n m 

- i) g(lK-Cj)tffc(k) = 2 ^ Wi(m; i) 7yfe (m, j, k) + ^ C H (m, n; i) <5;yfc( m , n, j, k) , 

1 m m.n 

(5.5) 

and 
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It is obvious that under the transformation of 

Aj(iJ), 7yfc(iJ,k), <$ij W (i,j,k,l), q(m)wi(i), q(l)wj(j)w k (k), g(l)g(p)J 

-»• ja/%(i,j), a7y fc (i,j,k), aS ijk i(i, j, k, 1) , ag(m)^(i), ag(l)«5 3 -(j)«5 fc (k), ag(l)g(p] 

(5.7) 

where a is any positive dimensionless constant, ( 15. 2 j) through ( 15. 6p are invariant, along with all 
the constraints of equality and inequality discussed. This property of scaling invariance has cer- 
tain implications for homogeneous turbulence modeling on the basis of objective optimization 
like It or Id] Specifically, it rules out the possibility that, under a given Vij ^ 0, a nontriv- 
ial bounded steady state solution is independent of the initial condition for the correlations. 
Otherwise, say, {/^ (i, j) , 7/^(1, j, k), 6lj kl (i, j, k, 1)} were such a solution, then for any a > 1, 
{a (i, j) , a7^ fc (i, j, k), a5^ kl {i y j, k, 1)} would be a solution too, and thus, the optimization 
scheme would produce an unbounded solution with the unboundedly great It and Id- 

We may infer from the above argument that solutions of the present homogeneous turbulence 
model for 



Ai(iJ), 7jfc(iJ,k), <% feZ (i,j,k,l) , g(m)tUi(i), g(l)w j (j)w A: (k), g(l)g(p] 

are time-dependent and initial condition dependent. This time-dependence appears compatible 
with the asymptotic behavior of homogeneous shear turbulence which shows the exponential 
growth of WkWk and so on, from the data of DNS and experiments (cf. [ID]). To explore this 
issue further, we take 

= S5 a 5 j2 , S>0 (5.8) 
and introduce an asymptotic form solution of 
Aj(iJ) = Pif{i,j)exp(aSt), 7 jfc (i,j,k) = S 7<g(i, j, k) exp(aSt), 

<5y«(i,j,k,l) = S' 2 ^(i,j,k,l)exp(a5't), g(m)u5i(i) = S'g(m)u) i (i) (0) exp(aS't), 



g(m)^(i) ^(j) = g(m)^(i) ^(j) (0) exp(a^), q(\)q(p) = S 2 g(l)g(p) (0) exp(a^) (5.9) 

where the quantities with superscript "(0)" are independent of time and a is a constant yet to 
be fixed. To see whether such a solution is admissible in the present model, we first substitute 
(Ojl and O into through (EH]) to get 

^(U) - ^E^^SVj) + E +TF(m;i)5a)A5 )) (m,j) 

m m 

- ^E D ( m ;M ( °W) + E (^(m,j) +F(m;j)« jl )j5f (m, i) 

m m 

= - X^/(m; n, i)7,^(n,m,j) - ^^(m; i) ^(m)^^)^ 

m.n m 

- J^A i (m;n,j)7,5J(n,m,i) - ^(m; j) g(m)wi(i) (0) , (5.10) 
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-^D(m;p)g(m)«5i(i) (0) = 2 W fc (m; p) /^(m, i) + ^ <^(m, n; p) 7j g](m, n, i) , 

m m m,n 

(5.11) 



<r®i J, k) - | J] D ^ *) 7$(m, j, k) + £ (V 4 ,(m, i) + W(m; i) ^) 7 $(m, j, k) 



m m 

r(0) 



- 5 E D ( m ^j) %-2( m > k > + E (^(mJ) + ^( m ; j) ^) 7 ^ (m ' k> i} 

m m 

- | E D ( m 5 k ) T*2( m > i. j) + E (^«( m > k ) + k) ^) 7 ^(m, i, j) 



m,n m 

A t (m; n, j) ^(n, m, k, i) - ^ B j( m '^) g(m)u) fc (k)«) i (i) ( ° } 

m,n m 

A z (m; n, k) ^(n, m, i, j) ~ E B k( m ' k ) g(m)«Ji(i) , (5.12) 



- D & i) g(lK-(j)^(k) (0) = 2 £ W,(m; i) 7 $(m, j, k) + ^ ^(m, n; i) 4%(m, n, j, k) , 

1 m m.n 

(5.13) 

and 

- E D & *) ^( 1 )^(p) (0) = 2 E ^( m ; i) g(p)^(m) (0) + C ji( m > n ; g(p)t2f( m )™i( n ) ° • 

1 m m,n 

(5.14) 

Here 

W(m; k) := i W(m; k) , W^m; k) := ^ Wi(m; k) , Wy (m, k) := i W^-(m, k) (5.15) 

which are independent of 5*. As expected from the linear structure of (15. 2 p through (15. 6p . 
the algebraic equations above are independent of time. Next, we can verify directly that, 
under ( 15. 9 p , all the constraints of equality and inequality listed in Section [2] and in (15. ip can 
be converted into the corresponding constraints of the same form on /3^' (i,j), 7 }^(i,j,k), 
^ijliih j) k, 1), etc. without the presence of S. 

The values of a, 0^ (i,j), 7 ^(i,j,k), <5^(i, j, k, 1) and so on, depend on v/S and L. To 
determine these values, we resort to maximize, say It = ^ k /3fcfc(k, k); Due to the supposed 
constancy of a and the directional sequence of t, we will maximize 

It = E^( k ' k ) | t=0 = E^( k ' k ) (5-16) 

k k 
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under (I5.10P through (I5.14p and the resultant constraints of equality and inequality mentioned 
above. Next, within the peculiar context of the asymptotic form solution of (15. 9p which holds 
mathematically for all t G (—00, +00), the flexibility of choosing any instant of time as the 
initial instant and the corresponding value change of the correlations obeying (15.91) motivates 
us to fix a specific upper bound for It above such as 

$><°>(k,k)</<°> or 5>£>(k,k) = /f. (5.17) 

k k 

Further detailed analysis and computation will check whether the solution matches with the 
data of DNS and experiments. Also, one can see the importance of developing an optimization 
algorithm and code so as to obtain the solution. 

One important issue yet to be resolved is whether, under a set of arbitrary initial conditions, 
a solution of (15. 2 p through (I5.6P exhibits the asymptotic behavior of (15. 9p . without the explicit 
constraint of the latter. 

As an apparently simple extension, one may also consider similarly the asymptotic solution 
under Vij = S(8n5j2 + <5i2<5ji) or Vij = S^SnSji — ^2^2) • For more general homogeneous turbu- 
lence of V^, instead of (15. 9p . we may seek an asymptotic solution of the form <p = (p^ exp(crt). 

6 Remarks 

We have discussed averaged turbulence modeling of multi-spatial scales for a turbulent flow 
of an incompressible Newtonian fluid in a bounded domain; the major points are summarized 
below. 

i. The Reynolds average is supposedly adequate to partition the flow fields into the aver- 
age ones and the corresponding turbulent fluctuations. There exists a function basis to 
decompose the fluctuations into the components associated with various spatial scales. 
Supposedly there is a probability density function / of these fluctuation components 
and there are well-defined statistical correlations of multi-scales among these fluctuation 
components (up to the fourth order). 

ii. The constraints of equality and inequality are formulated for the correlations in the wave 
number and physical spaces, which are either general intrinsic to turbulence or specific 
for the flow; the Cauchy-Schwarz inequality plays a significant role here. 

iii. The maximum information principle is resorted to determine the unbiased form of / un- 
der the equations of motion and the constraints of equality and inequality. Considering 
the formidable difficulty to implement the maximization of the information both ana- 
lytical and computational, closely related but simple alternatives like the invariants of 
the covariance matrix j)] are sought, two single objectives of special interest are, 
respectively, the determinant and the trace of the covariance matrix, the former suitable 
to the flow being slightly non-Gaussian and the latter due to its relation to the total 
fluctuation energy possessed in the turbulent flow. 

iv. The model with the multi-scale correlations up to the fourth order is the smallest model 
capable of resolving the pressure fluctuation correlation, which does not contain additional 
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coefficients besides the single objective to be maximized. Some preliminary results of its 
application to homogeneous turbulence are presented and the result on homogeneous 
shear turbulence suggests its possible compatibility with DNS and experimental data. 

To the goal of making the present formulation applicable, there are still many basic issues 
to be resolved, some of which are sketched below. 

1. It needs to be tested whether the information-based objective optimization reflects the 
physically observed phenomena. To check whether and how the model can produce mean- 
ingful results, we need to apply it to specific problems. Due to the number of variables 
and equations involved, it is preferable to consider some relatively simple problems first, 
like those of 1-D Burgers' equation and 2-D and 3-D homogeneous turbulence, with the 
help of various function bases. In this regard, the problem of homogeneous turbulence 
formulated in Section [5] will be investigated further. 

2. From the viewpoint of computation and analysis, a reduced model containing the correla- 
tions only of the second order (and at most the third order) has a great advantage, which 
is also an extension of the Reynolds stress equation model of engineering turbulence. The 
task is then how to model adequately the third order correlation related terms on the 
right-hand side of fl2.12j) . in the form of inequalities involving /^-(ij) and 7jj fc (i,j,k), so 
as to make the reduced model determinate. [12] attempted to address this issue, but 
the special tentative expression proposed is inadequate, for example, to satisfy the many 
constraints of equality and inequality involving the second order correlations listed in this 
work. 

3. The treatments of representative wave numbers, locality, and local isotropy at high wave 
numbers are expected to reduce the number of correlations and equations involved in a 
model. Moreover, these treatments are motivated physically, and due to their physical 
nature, they supposedly help narrow the spread-out of / and make the predictions better. 
How to formulate the treatments is yet to be resolved, especially in the need to satisfy 
the constraints on the correlations, particularly those holding in the physical space. 

4. Some simple approximate single objective functions are presented, alternative to the in- 
formation. Which approximate objective should be adopted is to be studied with the help 
of grid turbulence, etc. 

5. Pipe and channel flows of simple geometry will be simulated in order to test whether the 
model is adequate for wall-bounded turbulence. 

6. An algorithm of optimization and solution is needed and essential which can deal with 
the great number of constraints. 
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